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In this note, I summarise and comment on joint work with C. Bernardin, V. Kan- 
nan, and J. L. Lebowitz concerning two harmonic systems with bulk noises whose 
nonequilibrium steady states (NESS) are nearly identical (they share the same ther- 
mal conductivity and two-point function), but whose hydrodynamic properties (con- 
qq : vergence towards the NESS) are very different. The goal is to discuss the results in 

\Q | the general context of nonequilibrium properties of dynamical systems, in particular, 

what they tell us about possible effective models, or predictive approximations, for 
such systems. 
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1 Introduction 



This note is mainly based on the results derived in [Q] and summarised from a more physi- 
cal point of view in [2J. Both are made in collaboration with C. Bernardin, V. Kannan, and 
J. L. Lebowitz; hence the discussion here relies heavily on the work of others, although 
naturally the responsibility for any possible faults is solely mine. 

Two dynamical systems are considered in these works, combining the same harmonic 
Hamiltonian evolution with two different types of bulk noise. We consider d-dimensional 
crystals with the harmonic evolution determined by nearest neighbour interactions on a 
square lattice with fixed boundary conditions in the first direction and periodic boundary 
conditions in all other directions (if d > 2). The geometry is chosen mainly for compu- 
tational convenience. It has two natural boundaries, called the two ends of the crystal, 
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formed out of the sites attached to the fixed boundary particles. We denote the number of 
particles in the first direction by L and use the lattice spacing as the microscopic length 
unit: for instance, the microscopic volume is equal to the number of particles, iV = L d . 
By a thermodynamic limit of such systems, I mean taking L — > oo. 

Both of the above models have been studied before: the first one uses the self-consistent 
heat baths as in OlHl while the second velocity flip model has been considered in [|51|6ll21- 
In the self-consistent model, each particle is coupled to a Langevin heat bath and the in- 
put temperatures of the baths can vary along the chain but are uniquely determined by 
fixing the temperatures on the left and right end of the system at some given values, Tl 
and Tr, and then requiring that in the remaining sites there is no energy flux on average 
between a particle and its heat bath in the stationary state. As shown in [4], this require- 
ment leads to a system with a unique nonequilibrium steady state (NESS) satisfying the 
stationary Fourier's law. We also obtain explicit formulae for the heat conductivity and 
for the dependence of correlations on the self-consistent temperature profile. 

The bulk noise in the self-consistent model only conserves energy on average, and 
even then only in the NESS. The second bulk noise acts less disruptively, by randomly 
flipping the signs of particle velocities, and thus conserving energy in each realisation 
of the noise. More precisely, one can imagine that each particle carries its own (Poisso- 
nian) clock whose rings will enforce a velocity reversal, independently of what the other 
particles are doing. We then attach Langevin heat baths at the ends of the crystal, just 
as in the self-consistent system above. It is shown in Q that this model converges to a 
unique NESS which has a covariance matrix and thermal conductivity identical to that 
of the self-consistent model, as long as the parameters of the harmonic evolution and the 
boundary conditions match. 

Therefore, it might sound reasonable to conjecture that the thermal conduction prop- 
erties of the two systems are identical, even though both noises are "strong perturbations" 
in the sense that they convert the pure harmonic system with infinite conductivity [8] into 
a normally conducting one. In fact, our results strongly indicate that the NESS:s of these 
two systems are locally indistinguishable from each other: the difference between the 
NESS expectations of any local observable vanishes in the thermodynamic limit. How- 
ever, they are not globally identical since, for instance, the total energy fluctuations in 
the NESS differ by an amount which diverges as L — > oo. This turns out to be due to 
long-range corrections to local thermal equilibrium expectation values which are found in 
the velocity flip model but which are absent from the self-consistent model. Even more 
pronounced is the difference in the convergence towards the NESS: for the self-consistent 
model with pinning this occurs microscopically fast (exponential decay which is 0(1) 
in the original microscopic time scale) whereas it will typically take 0(N 2 ) microscopic 
time units before the state of the velocity flip model is close to its NESS. 

The purpose of this note is to recall the arguments why the latter behaviour should 
be the one most commonly found in physical systems, and to discuss how our earlier 
results can be interpreted as effective models for the energy transport in these systems. Its 
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outline follows closely a talk given by the author at a workshop in the NORDITA program 
Foundations and Applications of Non-equilibrium Statistical Mechanics in 201 1. Here I 
will try to clarify the key concepts and terminology in two explanatory sections, |2] and HI 
preceding the related results. 

2 On local thermal equilibrium 

The intuitive meaning of the terms local and global is clear: "local" refers to the prop- 
erties of the system in a (microscopic) neighbourhood of a point and "global" concerns 
the properties of the entire system. The distinction becomes relevant only for large sys- 
tems, when the entire system spans a large spatial region and these two length scales thus 
separate; in other words, in the thermodynamic limit. 

A local equilibrium state should thus be a state in which the properties of the system 
near a given point x in space are determined by some equilibrium state. For the present 
purpose, I will only consider equilibrium states determined by "canonical Gibbs mea- 
sures": given a collection of n conserved quantities Oi(q,p), i — 1, . . . , n, q,p £ (M. d ) N , 
the corresponding canonical Gibbs measures are labelled by n real numbers \ such that 
the measure e^ iXi0i( - q ' p Mqdp is normalisable. As long as the Lebesgue measure dqdp is 
invariant under the time-evolution — which is the case for any Hamiltonian evolution, also 
together with the above velocity flips — any of such Gibbs measures is stationary. Ther- 
mal equilibrium states are then defined by choosing as the observables the total energy, 
Oi(q,p) '■= H(q,p), and any other thermodynamically relevant conserved quantities, 
such as particle number, total momentum, etc. Temperature is then defined by T := (3~ l 
where (3 :— —X 1 typically needs to be positive for the measure to be normalisable. These 
states also serve as good reference measures of the self-consistent heat bath model in the 
sense that if all baths have the same temperature then the corresponding Gibbs measure 
is invariant under the time-evolution. 

We say that a state /i of the system is in local /z -equilibrium at x, if /i is an equi- 
librium state of the (infinite) system and (F(q,p)} )1 w (F(q,p)) IMj for all observables F 
localised at x. More precisely, this should hold for any F which depends only on qj,Pj 
for those j with \x — Xj\ < R where R > is a large microscopic length and Xj is the 
spatial location of particle j (usually, Xj = qj but for our crystals we use Xj = j). The 
correction term should vanish in the thermodynamic limit, and R can then be arbitrary as 
long as it is L-independent. The state ji is a local equilibrium state if this property holds 
at every point x £ R d . 

Thus to have a local equilibrium state necessarily implies the existence of a collection 
Po(x) of equilibrium states labelled by points x, although in the present generality, fxo(x) 
are not uniquely determined. If the reference equilibrium states are chosen from the 
thermal equilibrium states, then the state is said to be in local thermal equilibrium (LTE) 
and the distribution of the parameters Xi(x) is typically unique up to errors which vanish 
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in the thermodynamic limit. In particular, if we allow some thermodynamically small 
errors to the temperature distribution T(x), its meaning is then uniquely determined for 
any LTE system. Note that although then always T(j) = (p|) for any LTE state of a 
standard Hamiltonian system for which energy is the only conserved quantity, the above 
definition implies more: then necessarily also T(j + r) rj T(j) for any microscopic 
lattice displacement r. Thus LTE implies that T(x) must be essentially constant in any 
microscopic region of such systems. 

LTE states are common in simulations and in strongly stochastic particle systems, 
and it even seems reasonable to conjecture that in great generality an initial state will 
converge into some LTE state in a finite microscopic time. Indeed, it is questionable 
if the term "temperature distribution" should ever be applied to a state which is not at 
LTE. However, rigorous proofs of such generic convergence towards LTE have only been 
achieved in certain special cases, and it could well be the hardest part of a mathematical 
proof of dynamical Fourier's law. For instance, we have a proof of LTE in the present 
self-consistent model, which has a harmonic Hamiltonian term, but it remains the only 
part missing from an analogous proof when the Hamiltonian is anharmonic [|9). 

3 Local and global correlations: how can something be 
both vanishing and essential? 

How do the correlations behave in the two models described in the introduction? For 
simplicity, I will only consider the one-dimensional case in detail; some of these results 
immediately generalise to higher dimensions, as explained in [Q]|. Thus, from now on, 
L — N and d—1. 

The Hamiltonian part of the dynamics is harmonic, and thus by itself it has order iV 
conserved quantities and it transports energy ballistically through the system. In particu- 
lar, it has an infinite thermal conductivity. Our two bulk noises have been chosen precisely 
so that they would break most of these conservation laws without breaking the conserva- 
tion of the total (Hamiltonian) energy of the system. In the self-consistent case the energy 
is conserved only on average in the NESS (this is achieved by the self-consistent tuning 
of the bulk heat bath temperatures) and in the velocity flip model it is conserved in the 
bulk with probability one. However, as will be seen shortly, also other relevant conserved 
quantities can appear. 

As mentioned in the introduction, both systems have a unique NESS, and the two- 
point functions of these states, i.e., the covariance matrices of the variables {qj,pj}, coin- 
cide. The NESS of the self-consistent model is Gaussian, and thus uniquely determined 
by the covariance matrix. The NESS of the velocity flip model has some nonzero fourth 
order cumulants and thus cannot be Gaussian. Nevertheless, our results also indicate that 
these higher order cumulants all approach zero as iV — >■ oo. If this is true, all local corre- 
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lations would agree in the thermodynamic limit, and it follows from the known properties 
of the self-consistent case that both NESS:s must then be LTE states with the same tem- 
perature distribution. 

However, the convergence towards equilibrium is radically different in these two mod- 
els. In the self-consistent case with pinning, whatever the initial state, all local expecta- 
tions relax to their NESS value exponentially fast on a microscopic time-scale, while in 
the velocity flip case local energy needs a time 0(N 2 ) to equilibrate. In addition, even 
though they approach zero in the thermodynamic limit, the above mentioned corrections 
to LTE do contribute towards energy fluctuations in the NESS (this can happen since the 
energy fluctuations are determined by a sum over 0(N 2 ) local correlation functions). The 
next subsection summarises our findings. 

3.1 Comparison of energy fluctuations at the NESS 

In an infinite system without boundary terms, the bulk dynamics of the velocity flip model 
conserves the energy density, H(p, q; N)/N, N — > oo. The Hamiltonian function of the 
finite chain is given by 



where £j = £j(p, q) denote the local energy observables. Here we say that the system has 
pinning whenever v ^ 0. Without pinning (y = 0), also the total deformation density, 
N^ 1 r j> r j '■= Qj+i ~ Qj> is conserved and can fluctuate. (Since it is always, also with 
pinning, equal to q^+i — qo, it is trivially conserved in the finite system by the chosen 
boundary conditions. However, the point is whether its value can vary in the infinite 
volume canonical Gibbs state: this is possible if v = but not if v ^ 0.) 

In contrast, the bulk dynamics of the self-consistent model does not preserve the en- 
ergy density, although without pinning it does preserve the deformation density. This 
model is explicit enough to be studied rigorously, even in the thermodynamic limit: I 
quote below some of the results proven in [1] (without pinning) and in [@]| (with pinning). 

The correlations at NESS depend linearly on the input temperatures of the heat baths, 
(XiYj) = Yln=i Bxy(h j)T n , X, Y G {q,p}- Here B decays for increasing M : = 
1 + \i ~ j\ + \i — n\ + \j — n\, at least as fast as dictated by the following bounds: 

With pinning, (i, j) \ < ce~ aAI for some c, a > independent of N. 
Without pinning, for some iV-independent c > and denoting rj = qj + i — qj, 




j=0 



i:\i-j\ = l 



B%\i,j)\ < cM- 2 (l + lnM), \B%\i,j)\ < cM" 3 (l + lnM), 
S^(z,j)|<cM- 4 . 
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These results also yield the following estimate on energy fluctuations in the self-consistent 
NESS, 

1 1 N+1 

N< H '^ = N E^^) (eql ' Tj) + ^' 
j,f=o 

where we denote (A; B) = Cov(A, B) = (AB) - (A)(B) and 

with pinning R N = 0(N~%), 
without pinning R N = 0(7V-4ln 2 N). 

(The bounds for R N are likely off by a power N~? due to our poor control over the de- 
pendence of the self-consistent temperature profile on N.) Thus we can conclude that the 
energy fluctuations in the self-consistent model are O(N) and only the LTE expectations 
contribute to the leading behaviour. 

At the moment, we do not have any mathematically rigorous bounds for the energy 
fluctuations at the NESS of the velocity flip model, but numerical simulations give clear 
evidence that, although of same magnitude in N, they do not agree with the above LTE 
result. For instance, if Tl = 1 and Tr = 8, the LTE contribution to the N — > oo limit 
of sn '■= N (H '; H) / (H) 2 can be rounded to 1.20, while simulations produce 1.40(2) 
consistently for both N = 200 and N = 400 (more results are given in [2]). A possible 
mechanism, yielding an approximate value 1.40 in this case, is explained in Section[5J 

4 On kinetic and hydrodynamic scaling limits 

Various scaling limits have turned out to be useful in the study of conduction properties 
of dynamical systems. They share the feature that by considering a restricted collection 
of observables and an "unphysical" limit for some of the parameters of the dynamics, 
one obtains more easily solvable evolution equations for the limit observables. Therefore, 
such a limit produces an "effective theory" whose appropriately rescaled solutions yield 
approximations to the chosen observables of the original system; this typically under 
circumstances where a direct solution of the original system is not possible. 

The term "effective theory" is commonly used in physics but, in my experience, it 
takes some effort to explain the concept to those unfamiliar with it. After all, it is "effec- 
tive", not because it is very accurate (usually rather the opposite) but because it bypasses 
some obstacle, and it is not really a theory at all but a model for some parts of the original 
"theory". Personally, I would prefer to replace this term with something more immedi- 
ately descriptive, such as predictive approximation: the result is an approximation which 
can be used to predict some otherwise uncontrollable properties of the original system. 

Examples of such scaling limits abound: mean field limits, adiabatic limits, renormal- 
isation group transformations, etc. Here I will only discuss two such limits in detail: the 
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hydrodynamic and kinetic scaling limits. In a hydrodynamic scaling limit space is scaled 
by and time by L 2 . For nonequilibrium problems L is usually chosen so that the size 
of the system remains finite in the limit, and thus it can have boundaries through which 
the nonequilibrium state can be generated. For instance, a solution o(x', t') of a predictive 
approximation could be related to a solution of the original system by 



for all large enough L. Here the observable Oj(t) := Oj(q(t),p(t)) measures some mi- 
croscopic property carried by the particle j, F is an arbitrary compactly supported test- 
function and Xj(t) := Xj(q(t),p(t)) denotes the spatial position of the particle j at time 
t. The hydrodynamic scaling limits are best suited to study diffusive phenomena, as the 
scaling leaves such invariant, and the resulting predictive approximations typically (and 
also here) involve diffusion processes. 

The hydrodynamic scaling limits I would like to contrast with kinetic scaling limits. 
Although the latter term is not as firmly established as the first one, it can be motived by 
two properties shared by the limit evolution equations: First, the scaling leaves velocities 
invariant and the evolution will typically be dominated by constant velocity, i.e., ballistic 
motion. Secondly, the resulting limit equations are often those found in kinetic theory, 
such as Boltzmann transport equations. 

A common way to arrive at a kinetic scaling limit is to start with a system which 
has some explicitly solvable "free" dynamics involving motion with constant velocity. 
Free classical particles form obviously such a system but so do many wave-equations, at 
least if one considers the evolution of the Wigner functions of the wave fields. Then a 
perturbation is introduced into the system, for instance, by adding a potential term XV. In 
any case, let A > denote the "strength" of the perturbation with A = corresponding 
to pure free evolution. The scaling limit is then defined by scaling space and time with 
A a , a > 0, and considering A — > 0; this obviously leaves all velocities invariant. Here 
different choices of a > typically produce different types of limit evolution. However, 
most often there is a unique choice a such that the limit equations are ballistic for any 
a with < a < a but become nonballistic for a = a . For instance, for many weakly 
perturbed wave equations a = 2 (for rigorous results, see e.g. lfT0l[TTl[T2l '). 

If the limit evolution is not diffusive, it is not clear how to use hydrodynamic scaling 
limits. At the very least one needs to modify the scaling functions, but how? Here ki- 
netic scaling limits can be very helpful by identifying the first nontrivial effects produced 
by the perturbation (for any time scale 0(A~ a ), a < a , the motion is ballistic, hence 
trivial). Then one can study what happens to the predictive approximations — produced 
by rescaling the solutions of the kinetic limit equation — at times longer than the kinetic 
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time scale. If no new effects appear, the overall result can be a simple correction to the 
constants predicted by the kinetic equation. This seems to be common for normally con- 
ducting systems where the diffusion constants computed from the hydrodynamic scaling 
limit of the kinetic equation determine the leading behaviour of conductivity as A — ^ 0. 
However, the terms neglected in the kinetic scaling limit can also alter the character of 
evolution entirely. For instance, the Boltzmann equation arising in the kinetic limit of 
the two-dimensional Anderson model is expected to be diffusive ITT3T1 while the original 
model should be exponentially localised, with the localisation length diverging for A — > 
(for related experimental evidence of such phenomena, see Ifl4l0 . 

Thus one should not treat the "effective theories" obtained from the scaling limits as 
totally universal, but rather as predictive approximations which might persist for times 
longer than those indicated by the scaling, but do not necessarily do so. One has to be 
particularly careful if the predictive approximation generates singularities. These singu- 
larities might conspire with some of the neglected terms and affect the evolution already 
during finite kinetic time scales. This appears to happen in Bose condensation in a bosonic 
Boltzmann-Nordheim equation [TT3TI . and the "entropy solutions" to Euler equation pro- 
vide another example of a nontrivial continuation beyond a singularity. 



5 Fluctuating hydrodynamics of the velocity flip model 

Only the velocity flip model without pinning will be discussed here since its hydrody- 
namic structure is richer than the one with pinning or in the self-consistent model. The 
local quantities building up the energy and deformation densities, as well as the associated 
currents, are then 

^,2 ^2 -i 

r j = q j+ i-q j , j) = ~Pj ■ 

Explicitly, the observables and their currents satisfy Ej (t) —Sj (0) = f* ds [j'J (s) — (s 

and r,(t) - r,-(0) = /> [j» - f j+1 (s)] . 

We parameterise the canonical Gibbs measures in terms of ft — T 1 > and rGl 



N+l 
3=0 



and denote expectations over this measure by (-)t,t- Then 

(p 2 J ) T ,r = T, (£j) TtT = T + r 2 /2, (r J ) T , T = r. (2) 

We next assume that LTE with respect to the above measures holds for all times. Then 
the mean energy density e(x', t') and deformation density u(x', £'), under hydrodynamic 
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scaling, satisfy e(x' } f) = T(Nx', N 2 t') + t{Nx', N 2 t') 2 /2 and u(x', f) = t(Nx', N 2 t'), 
up to small errors. For this particular model, there is also a happy accident, and the 
evolution equations of the hydrodynamic fields can be closed merely by assuming LTE: 
Since (£ denotes the generator of the stochastic process) 

f. = _(V0), + C(hj), j) = -T^Vr),-! + Cil^Pj), 

the slightly formal argument in [1] yields that for t > 0, x G [0, 1] the fields u(x, t) and 
e(x, t) should satisfy 

d t u = 1 - l d 2 x u, (3) 
d t e = (2 1 )- 1 dl(e + u 2 /2), (4) 

with the boundary conditions 

(d x u)(0,t) = (d x u)(l,t) = 0, (5) 
(e - u 2 /2) (0, t) = T L , (e - u 2 /2) (1, t) = T R . (6) 

In fact, in a recent preprint [|T6l an analogous result has been rigorously derived for 
slightly different boundary conditions and assuming that the initial state is close, in the 
sense of relative entropy, to an LTE state. 

Thus the above result already provides a predictive approximation for the diffusive 
scale {i.e., macroscopic) evolution of mean energy and deformation densities in this sys- 
tem. However, there are natural questions for which it is not directly predictive, the 
fluctuations of the total energy being one of them. To study the fluctuations of the hydro- 
dynamic fields around their expectation values e and u, we define, for an arbitrary choice 
of test functions F and G, the fluctuation field observables as 

u{j/N,t)] , 
e{j/N,t)] . 

This roughly coincides with the example given in Section 0] since here the particles are 
embedded at their lattice sites, i.e., we use Xj = j. The main differences are the miss- 
ing time-average and the different scaling of the fields. The above choice of scaling is 
determined by requiring that the fluctuation fields should have nontrivial limits. 

The fields TZ^ and yf can also be interpreted as time-dependent distributions. We 
argue in [|T) that TZf — > 1Z t and yf* — > y t , in the sense of distributions as N — > oo, and 



nf{F) = -=Y, F ^IN)HtN 2 )- 

V j = i 

1 - 

y t N (G) = -j=Y i GU/N) [S.itN 2 )- 

7=1 
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that the limit distributions solve the following stochastic differential equations: 

d t n = 1 - 1 d 2 x n-d x (cw 1 ) , (7) 



d t y = (27)- 1 {d 2 x (ull) + d 2 x y) - d x (cuW 1 + c^Tj2W 2 ) . (8) 

Here W\, W 2 are independent space-time white noises, T = T(x, t) = e(x, t) — \u(x, t) 2 , 
and c = c(x, t) = ^j2^~ 1 T(x, t). Combined together with the evolution equations for e 
and u, these results form another predictive approximation of the evolution of energy and 
deformation densities in the original system, however, with a greater "resolution" than 
if the fluctuation fields are neglected. These stochastic evolution equations are derived 
using similar techniques and assumptions as those for the hydrodynamic equations ©- 
©. Most notably, LTE is assumed to hold. 

Controlling the fluctuation fields makes it possible to predict also the total energy fluc- 
tuations at the NESS which can then be compared with the numerical simulation results 
mentioned in Section [3TTT As t — > 00, the solutions to the hydrodynamic equations sat- 
isfy u(x,t) — > and e(x,t) — > T(x) where T(x) denotes the linear profile connecting 
the boundary heat bath temperatures T L and T R . The fluctuation fields converge into two 
independent Gaussian fields, 1Z and y, with covariances 

(11(F) 2 ) = [ dxT(x)F 2 (x), 
Jo 

(y(G) 2 ) = f\xT(x) 2 G(x) 2 + (T L - T R ) 2 [ 1 dxG(x)((-A )- 1 G)(x) , 
Jo Jo 

where A denotes the Laplacian with Dirichlet boundary conditions on [0, 1]. (Somewhat 
surprisingly, an analogous computation with pinning yields the same NESS energy fluctu- 
ation field y.) We then find the following explicit prediction for Sqo, the thermodynamic 
limit of sn evaluated at the NESS, 

4T L r R + f(r L -r R ) 2 

(Tl + Tr)2 • w 

Inserting the boundary values used in the simulation thus yields Sqo ~ 1.403 in perfect 
agreement with the numerically observed value. This is probably our strongest evidence 
at the moment that the above fluctuating hydrodynamics indeed forms a predictive ap- 
proximation of the hydrodynamic properties of the original lattice system. 
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